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A GUIDE TO USING METEOROID-ENVIRONMENT MODELS 

FOR EXPERIMENT AND SPACECRAFT DESIGN APPLICATIONS 

By Donald J. Kessler 
Manned Spacecraft Center 

SUMMARY 


In this report, a method is derived for transforming a meteoroid flux (defined 
relative to the Earth and expressed as a function of meteoroid mass, velocity, and 
mass density) into a meteoroid flux expressed as a function of penetration thickness 
(or some other parameter, such as energy or momentum) on a moving spacecraft at 
some particular distance from the Earth. The meteoroid flux on a moving spacecraft 
is increased over the meteoroid flux relative to the Earth. Use of specific weighted- 
average velocities and mass densities is necessary for transforming the meteoroid 
flux- mass distribution relative to the Earth into a flux- penetration distribution on a 
moving spacecraft at a specified distance from the Earth. The procedure for perform- 
ing the transformation is summarized in this report, and an example is discussed. 


INTRODUCTION 


The conversion of a meteoroid distribution (expressed as a function of mass, 
velocity, and mass density) into a distribution expressed as a function of some other 
parameter can be accomplished by several methods. One method, which gives an 
approximate solution, is to average all meteoroid velocities and densities and convert 
the mass distribution into a distribution of the new parameter. In the use of such a 
method, the variables involved are not adequately considered. For example, if a 
small high-velocity meteoroid penetrates a material to the same thickness as a large 
low-velocity meteoroid, then flux as a function of penetration thickness could change 
as the velocity distribution changed, even though the average velocity might remain • 
constant. As is demonstrated in this report, the velocity distribution can be replaced 
by a weighted-average velocity, where the proper weighting is a function of the equa- 
tion describing the desired parameter, as well as a function of the slope of a log-flux- 
versus- log- mass curve. . 

The approximation is also usually made that the meteoroid environment relative 
to the Earth is the same as that for an orbiting spacecraft. Actually, both the relative 
meteoroid impact velocities and the meteoroid flux increase as the spacecraft velocity 
relative to the Earth increases. 



Methods are developed in this report to convert a meteoroid flux expressed as a 
function of mass, velocity, and mass density (relative to the Earth) to a flux expressed 
as a function of some other parameter relative to an orbiting spacecraft. It is assumed 
that the spacecraft has no preferred orientation with respect to the Sun and that the 
spacecraft orbit about the Earth has no preferred plane. It is shown in this report that 
this assumption is equivalent to assuming that the geocentric flux is not directional. 

The author wishes to thank Kenneth Baker, who first suggested an approach sim- 
ilar to the one taken in this paper; Herbert A. Zook, who assisted the author in nu- 
merous instances; and Burton G. Cour-Palais, who suggested that this work be organized 
and published. The author is also indebted to the other members of the Meteoroid 
Sciences Branch (dissolved October 1970) at the NASA Manned Spacecraft Center. 


SYMBOLS 


A normalization constant in the impact-angle distribution equation 

B constant coefficient in the penetration equation (or in the equation defining 
another parameter) 

C constant coefficient in the cumulative-flux-versus-mass equation 

D normalization constant in the equation for the meteoroid velocity distribution 
relative to a moving spacecraft 

d^ radius of an imaginary sphere concentric to and outside the Earth 

F cumulative meteoroid flux, number of impacts per unit area per unit time 

F cumulative meteoroid flux on a stationary spacecraft with respect to the Earth, 

e number of impacts per unit area per unit time 

F cumulative meteoroid flux on an orbiting spacecraft, number of impacts per unit 
area per unit time 

f ratio of the mass flux of meteoroids with velocity V m on an orbiting spacecraft 
to the mass flux of meteoroids with velocity V m on a hypothetical spacecraft 
at the same position and stationary relative to the Earth 

f, ratio of the mass flux of meteoroids with a velocity distribution of V on an 

t m 

orbiting spacecraft to the mass flux of meteoroids with a velocity distribution 
of V m on a hypothetical spacecraft at the same position and stationary with 

respect to the Earth 

G ratio of the meteoroid mass flux at a distance r from the Earth to the meteor- 
oid mass flux near the surface of the Earth (if Earth shielding is neglected and 
only meteoroid velocity is considered) 
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ratio of the meteoroid mass flux at a distance r from the Earth to the meteor- 
oid mass flux near the surface of the Earth (if Earth shielding is neglected and 
only a meteoroid velocity distribution of is considered) 

direction from which a meteoroid stream comes 
meteoroid mass 

power to which velocities are averaged 
number of meteoroids having the property ( ) 
distance from the center of the Earth 
radius of the Earth 

number density of meteoroids near the Earth 

thickness a meteoroid will penetrate (or exhibit some other behavior), depending 
on the values of j3, y, 6, and e 

impact velocity of a meteoroid on the spacecraft 

representative velocity of meteoroids when converting a flux expressed as a 
function of mass to a flux expressed as a function of t 

escape velocity from the Earth surface 


velocity of a meteoroid relative to the Earth before interaction with the Earth 
gravitational field 

relative velocity between the Earth and a meteoroid at some distance from the 
Earth 

velocity of an orbiting spacecraft relative to the Earth 


relative velocity between the Earth and a meteoroid at the edge of the Earth 
atmosphere 

exponent of m in the flux- versus- mass equation 
exponent of m in the penetration equation 
exponent of V in the penetration equation 
exponent of p in the penetration equation 
exponent of cos 9 in the penetration equation 


3 


tj Earth shielding factor resulting from ’’line of sight" Earth shielding 

r) effective Earth shielding factor resulting from the curved paths of meteoroids 

e under the influence of Earth gravity 

6 angle between the normal to a surface and the negative velocity vector of an 
impacting meteoroid 

6 representative impact angle of a meteoroid when converting a flux expressed as 
a function of mass to a flux expressed as a function of t 

A angle between the velocity vectors V m and V g 

p meteoroid mass density 

p representative meteoroid mass density when converting a flux expressed as a 
a function of mass to a flux expressed as a function of t 

0 angle between the normal to the surface of the Earth and the meteoroid negative 

velocity vector V„ 
m 

half angle of the Earth when viewed from a distance r 

4 * effective half angle of the Earth when viewed from a distance r, when the curved 
e trajectory of the meteoroid between r and the Earth is considered 


DIRECTIONALITY RELATIVE TO THE EARTH 


The meteoroid flux relative to the 
Earth is probably directional, although the 
degree of directionality is not yet well de- 
fined. However, if, from all points on the 
Earth, meteoroids were observed entering 
the Earth atmosphere, the total distribu- 
tion of meteoroid entry angles would be 
random, as is shown by the following 
analysis. 

A random distribution can be de- 
scribed mathematically as follows. Con- 
sider a unit-area flat plate sitting under a 
large hemispherical dome (fig. 1). Each 
element of surface area of the dome emits 
an equal number of meteoroids toward the 
flat plate. (That is, all meteoroid direc- 
tions are equally probable relative to the 
plate.) Let e be the angle between the 


i 



Figure 1. - Random flux on a flat plate. 
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normal to the flat plate and the direction from which the meteoroids are arriving. The 
number of meteoroids hitting the plate from angle 0 (±de /2) will be proportional to 
the element of dome surface area having angle 0 (i. e. , sin e ) and the component of 
plate surface area exposed to direction e (i.e., cos 0), or 


n (0 ) = A sin 0 cos 0 

0 


( 1 ) 


where n (0 ) is the number of meteoroids hitting the plate from angle 0 , and A is 

a normalization constant. If the total number of meteoroids hitting the plate is nor- 
malized to 1, then 


jr/2 

J (e)d0 = 1 (2) 


and the value of A is 2. 


Thus, the normalized random distribution of meteoroids is 


n (0) = 2 sin 9 cos e 

U 


(3) 


Zero-Gravity Earth 

Consider a group of meteoroids, of 
uniform flux throughout, which comes 
from direction i and hits the Earth. (This 
situation is similar to the Earth passing 
through a meteoroid stream, as in fig- 
ure 2. ) Neglect, for a moment, Earth 
gravity. If meteor- observation stations 
were located over the entire hemisphere 
exposed to direction i, the entry angles 
observed would be proportional to the 
area of the Earth where observation of an 
entry angle e was possible (i.e., sin 0), 
times the component of that area facing 
direction i (i.e., cos 0). Thus, the 
angular distribution observed by the total 
of all observers would be 



Meteoroid 
stream from 
direction i 


Figure 2. - Meteoroid stream hitting the 
Earth (gravity neglected). 


n . (0) = A. sin 0 cos 0 (4) 

0 j 1 1 
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If j directions (streams) existed, the total angular distribution observed over the en- 
tire Earth surface would be 


i=J 


i=l 


(5) 


or 


i=j 

n 0 ( 0 ) = sin 0 cos 0 A. 

i=l 


(6) 


If equation (6) is normalized so that the total number of meteoroids striking the Earth 
is 1, then 


i-j 


L 



(7) 


and 


n ( 0 ) = 2 sin 0 cos 0 (8) 

U 

Gravitating Earth 

Consider a meteoroid stream from direction i having a relative velocity V^, 

with respect to the Earth before interaction with the Earth gravitational field and a ve- 
locity upon impact with the Earth. Define an imaginary sphere of radius d^, con- 
centric to and outside of the Earth, such that, if the trajectory of a meteoroid in the 
stream were to pass the Earth at a distance greater than d^ (when gravity is ne- 
glected), then the gravitational force of the Earth would be insufficient to cause the 
meteoroid to hit the Earth. If the meteoroid trajectory (when gravity is neglected) 
were to carry the meteoroid a distance of exactly d G from the Earth, the gravity- 

affected meteoroid would just graze the Earth, and the meteoroid perigee would be the 
radius of the Earth r g (fig. 3). 


6 




Thus, from conservation of angular 
momentum (ref. 1) 


r e V ~ = d Q V G 


(9) 


Conservation of energy requires that all 
meteoroids with velocity will have 


Figure 3. - Meteoroid stream hitting the velocity V M when they strike the Earth. 
Barth (gravity Included). A trajectory (when gravity Is neglected) 

that passes the Earth within a distance d, 

will pass through the imaginary sphere of radius d^, and make an angle e between 

the normal to the surface of the imaginary sphere at the point of intersection and the * 
meteoroid trajectory. The meteoroid will hit the Earth at angle 0, which is measured 
between the normal to the Earth surface at the point of impact and the negative velocity. 


vector V . Conservation of angular momentum requires that 


r V. sin 0 = d_V 0 sin e 
e 00 G G 


( 10 ) 


If equations (9) and (10) are combined 


sin 0 = sin e 


( 11 ) 


From equation (4), the total distribution of the intersection angles between the 
meteoroid trajectories and the normals to the surface of the imaginary sphere at the 
respective points of intersection is 


n . (e)de = A. sin e cos e de (1*2) 

e, i i 


Thus, from equations (11) and (12), it becomes obvious that the total angular distribu- 
tion of the intersection angles between the meteoroid trajectories and the normals to 
the surface of the gravitating Earth at the respective points of impact is 


n, .(0)d 0 = A. sin 0 cos 0 d 0 
1 1 


(13) 
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If equation (13) is normalized 


n^(0) d0 = 2 sin <p cos 0 d0 


(14) 


Thus, the entry-angle distribution for meteoroids, averaged over the Earth sur- 
face, is random, regardless of the initial meteoroid directionality or the effects of 
Earth gravity. It does not follow, however, that the entry-angle distribution is random 
at some particular point on the Earth surface. 


TRANSFORMATION OF THE FLUX-VERSUS-MASS EQUATION 


If a meteoroid flux expressed as a function of meteoroid mass, velocity, impact 
angle, and mass density is given, the problem arises of converting these distributions 
to the meteoroid flux expressed as a function of some other parameter. Since the 
parameter usually required is penetration thickness, the remainder of the analysis 
will refer to penetration thickness as the parameter to which the transformation is to 
be made. However, the development will be general enough to include any other 
parameter, such as energy or momentum. 

Assume that the thickness t that a meteoroid will penetrate (over the range of 
possible values of m, V, p, and e ) can be expressed by 


t = Bm /3 V r p 6 cos € e 


(15) 


where m is the meteoroid mass, V is meteoroid velocity, p is meteoroid mass 
density, and e is the impact angle between the meteoroid trajectory and the normal 
to the spacecraft surface at the point of impact; B, 0, y, 6, and e are constants 
that depend on the structure of the spacecraft, the spacecraft penetration mode, and 
the particular parameter desired (such as penetration thickness). Within the limits of 
possible values of meteoroid velocity, mass density, and entrance angle, a range of 
meteoroid masses that will penetrate a thickness t exists. Assume that over this 
range of mass values, the mass distribution can be expressed as 


GL 1 

n (m)dm = aCm” " dm 


m 


(16) 


where n m (m) is the flux of meteoroids of mass m, and a and C are constants. 

The flux of meteoroids of mass m and larger is then given by integration of equa- 
tion (16) from m = m to m = 00 , or 


F = Cm’ ff 


(17) 
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It would be convenient to use a single representative meteoroid velocity V , a single 

Si 

representative meteoroid mass density p , and a single representative impact angle 

a 

e in equation (15) and to combine equation (15) with equation (17) to obtain 


F = CB a// ^t~ a ^v v a /Pp 5a// P C os ea ^e 

ci Si Si 


(18) 


The problem is thus reduced to defining V , p , and 8 so that equation (18) is 
correct. a a a 

Assume that the flux of meteoroids with mass m, velocity V, mass density p, 
and impact angle e can be expressed as 


n t (m,V,p, e ) = n m (m)n v (V)n p (p)n 0 ( 0 ) 


(19) 


where Oy(V), n ^(p), and n^ ( 0 ) are the meteoroid velocity distribution, mass-density 

distribution, and impact-angle distribution, respectively. Each distribution, except 
n m (m), is found by comparing meteoroids of the same mass and is normalized so that 


J n v (V)dV = 1 
v=o 


OO 

f n p (p)dp = 1 

P= 0 



n 0 (e)de = 1 


( 20 ) 


( 21 ) 


( 22 ) 


The flux that will penetrate to thickness t and greater is then (ref. 2) 


.77/2 


F= / n 9 ,„ / n p (p) / "yW / n m (m)dm dV dp de (23) 

p= 0 V=0 m=m. 


e=0 
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where m t is the mass of the meteoroid that will penetrate thickness t at velocity V, 

mass density p, and impact angle e. The value of m^. is found from equation (15) 
and is expressed as 


m. 


= ,%iV/V s/ ^ e 


(24) 


Thus, by using the expression for n (m) in equation (16) and integrating equation (23) 
over mass 


F = CB a/ V a//3 


f n^(g) f n p(p) f n y(V)V^ a ^p 6a//, ^cos ea ^0 dV dp de 

9=0 p =0 P v=o 


(25) 


Since the distributions of V, p, and e are normalized (eqs. (20), (21), and 
(22)), then by definition 


f n v (V)V ya///3 dV s Y ya ^ (26) 

v=o 

f n (p)p 6a/ ^dp - p 6a/ 13 (27) 

p=0 



n ( 9 ) cos ea ^e dg = cos ea ^9 
0 


(28) 


Thus, equation (25) becomes 


F = CB a/0 t'“ /S V W(3 p ia/fS cos £ “ /f! 9 


(29) 
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If equations (18) and (29) are compared 


v a =G^ /yc 

Pa ^r 


cos e 

a 



(30) 

(31) 

(32) 


Thus, the representative values to be used in equation (18) are defined in equations (30), 
(31), and (32). 

The environments defined in references 3 and 4 can be used as examples for 
using equations (30), (31), and (32). As developed in the section entitled Gravitating 
Earth, the impact-angle distribution (relative to the entire Earth surface) is random. 
Thus, for a random directional flux (eq. (8)), equation (32) becomes 


ea/i 3 

cos ' r e ~ 
a 


r ^/2 

= / 2 sin e cos e 

e=0 


ea/p 
cos 'e 


de 


(33) 


or 


cos 0 

a 



(34) 
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Power to which velocities are averaged, n 


Figure 4. - Weighted-average meteoroid 
velocity relative to the Earth as a func- 
tion of weighting (velocity distribution 
from ref. 3). 


The velocity distribution given in 
reference 3 was used to compute the 

weighted-average velocity \V / for 
various values of n. The results are 
shown in figure 4. If spacecraft shielding 
were being designed relative to this veloc- 
ity distribution, then 



(35) 


However, as is developed in the section 
entitled Meteoroid Flux and Impact- 
Velocity Distribution for an Orbiting 
Spacecraft, such a velocity distribution is 
modified relative to an orbiting spacecraft. 

Since a meteoroid mass-density dis- 
tribution is not given in reference 3, then 
the average mass density given in refer- 
ence 3 is also the representative average, 
or 



However, in the mass-density distribution given in reference 4, 51 percent of the me- 

3 

teoroids have a mass density of 0. 37 g/cm , 45 percent have a mass density of 

3 3 

2. 8 g/cm , and 4 percent have a mass density of 8 g/cm . Thus, from this mass- 

density distribution 


p a = (o. 51 x 0.37 6a///3 + 0.45 x 2. 8 6a//3 + 0. 04 x 8 ba ^Y' ^g/cm 3 (37) 


METEOROID FLUX AND IMPACT- VELOCITY DISTRIBUTION 
FOR AN ORBITING SPACECRAFT 

Thus far in spacecraft design, the assumption has been ms •- lhat either the me- 
teoroid hazard for an orbiting spacecraft is the same as the meteoroid hazard relative 
to the Earth, or the meteoroid environment is already defined relative to an orbiting 
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spacecraft. However, this is rarely the case. Meteoroid fluxes and velocity distribu- 
tions obtained from meteor observations are, by the nature of the observations, rela- 
tive to the Earth. Although meteoroid fluxes measured by Earth satellites (ref. 4) are 
relative to an orbiting spacecraft, the velocity distribution used to reduce the 
penetration- thickness data is obtained from photographic meteors. Thus, mass deter- 
minations from the penetration thicknesses are made from a velocity distribution rela- 
tive to the Earth. Such approximations do not introduce large errors, but the nature of 
these errors should be known. In this section, the problem of transforming the mete- 
oroid flux and velocity distribution relative to the Earth to a flux and velocity distribu- 
tion relative to an orbiting spacecraft is discussed. 


Velocity Distribution From Uniform-Velocity Meteoroids 

Consider a spacecraft in orbit about the Earth with a velocity vector V g rela- 
tive to the Earth. Consider also a group of meteoroids, each with a velocity vector 

V relative to the Earth. The velocity of any one of these meteoroids relative to the 
m 

spacecraft is given by the vector relationship 


which can also be expressed as 


V = V - V 
m s 


V 2 = V 2 + V 2 - 2V V cos X 
m s ms 


(38) 


(39) 


where X is the angle between the velocity 

vectors V and V , as shown in figure 5. 
ms’ & 

Thus, a meteoroid velocity distri- 
bution develops relative to the spacecraft 
(depending on the angle X) from a single 
meteoroid velocity interval relative to the 
Earth. As is shown in the section entitled 
Gravitating Earth, if the spacecraft 
spends equal time over all parts of the 

Earth, all directions for V m are equally 

probable, the exception being the direc- 
tion shielded by the Earth. However, for 
the purpose of this analysis, Earth shield- 
ing will be ignored. 



Figure 5. - Vector relationship between 
the spacecraft velocity relative to the 
Earth and the meteoroid velocity rela- 
tive to the Earth. 
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Thus, the number of meteoroids coming from direction A will be proportional 

to the number of directions V will point as it rotates around the V -axis or 

m s 


n^(A) rc sin A 


(40) 


Since the probability that the spacecraft will encounter meteoroids increases with V, 
the meteoroid velocity distribution relative to the spacecraft is also proportional to V. 
(For example, if V = 0, the spacecraft never encounters a meteoroid. ) Therefore 


n g (V)dV °c V sin A dA 


(41) 


If equation (39) is differentiated 


sin A dA = 


V dV 

V V 
m s 


(42) 


Thus, if equations (41) and (42) are combined 


n g (V)dV = DV 2 dV 


(43) 


where V is between V - V and V + V , and D is the normalization constant. 

ms ms 

If the assumption is made that V < V , the value of D is evaluated by normal- 
izing the velocity distribution so that m 


Thus 



+V 

m s 



n g (V)dV = 1 



(44) 


(45) 
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and 


n (V)dV = 
s 


3V 2 dV 

2V (3V 2 +V 2 ) 
s\ m s / 


( 46 ) 


where V -V s V< V +V, and V s V . The function n (V)dV is then the 
m s m s s m s 

fraction of the meteoroids which impacts the spacecraft with a velocity between V and 

V + dV and which results from meteoroids with velocity V m relative to the Earth. 

If V sV , then by exchanging signs on the lower limit of integration in equa- 
tion (44) m s 


_ _ 3V 2 dV 

n (V)dV - 7 -~2 

2V m (3V c Z + V m ) 

mV s m / 


(47) 


where V - V <V^V +V, and V ^ V . The weighted-average velocity 
sm ms’ ms 


fc= v 


1 /n 


/ is found from equation (46) (assuming that in most cases V s V ). 

S III 


/ 


v +v o 

m s 


— V -V 
^.n _ m s 


V n n (V)dV 
s 


/ 


v +v 

m s 


(48) 


V -V 
m s 


n g (V)dV 


or 



(49) 
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where V g < V m> and the functional dependence is introduced to remind the reader of 

the V dependence, 
m 

An example of a direct application of equation (49) follows from considering a 
spacecraft in orbit about one of the large planets (i. e. , Jupiter, Saturn, Uranus, or 
Neptune). The escape velocities of these planets are so large (compared with their 
orbital velocities around the Sun) that all meteoroid velocities near the planet are es- 
sentially the escape velocity of that planet. Thus, if the escape velocity of one of the 
planets is unity (i. e. , V m = 1), and if a spacecraft is in circular orbit around the 

planet, V g = 1//2. If n = 1, equation (49) gives an average impact velocity of 1. 29, 

which is 29 percent higher than the meteoroid velocity relative to the planet. If n = 2, 
equation (49) gives a root- mean- square velocity of 1. 33. 

Before equation (49) can be used with a distribution of V , a flux increase fac- 
in' 

tor, which is also a function of V and V , must be determined. This factor will be 
derived in the following section. s m 


Increase in Meteoroid Flux on an Earth-Orbiting Spacecraft 

The increase in meteoroid velocity relative to an orbiting spacecraft (discussed 
in the previous section) also leads to an increase in the meteoroid flux on the orbiting 
spacecraft, as compared with the meteoroid flux and velocity detection on a hypotheti- 
cal spacecraft located at the position of the orbiting spacecraft, but which is stationary 
with respect to the planet. If the number density of meteoroids near the Earth is S, 
the meteoroid flux detected on a spacecraft near the Earth and which is stationary with 
respect to the Earth is given in reference 5 as 


F 


_ 1 
e 4 



(50) 


where V m is the meteoroid velocity relative to the Earth. Similarly, the meteoroid 
flux on an orbiting spacecraft will be 


F 


_ 1 
s 4 



(51) 


where V is the meteoroid velocity relative to the orbiting spacecraft. 
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If only one group of meteoroids, those with velocity V , is considered, then the 

ratio of the meteoroid flux on the orbiting spacecraft to the flux on the stationary space- 
craft becomes 



where n (V) is given by equation (46) or (47), and the limits of integration are 
s 

V - V and V + V for equation (46) or V - V and V + V for equation (47). 
ms ms n s m ms v/ 

The flux increase factor then becomes 



3V 2 + V 2 
m s 


3V 


for V < V 
s m 


m 


(53) 


and 


3V 2 + V 2 

<V m ) = - for V_ > 


3V V 
s m 


s m 


(54) 


The total flux increase, resulting from all values of V , is then given by 


>t =/ fn V, m( V m) dV m 


(55) 


where m. (V ) is a fraction of meteoroids having velocity V . 

v , m' m m 

If V < V , which is the case for a spacecraft in orbit around the Earth, then 
s m 


f t =1 + 


v 2 v " 2 

s m 
3 


(56) 
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For example, figure 4 shows that V 


-2 


m 


- 1/2 


has a value of 16. 8 km/sec for the 


velocity distribution given in reference 3. A spacecraft in circular orbit just above 

the Earth will have a V value of 7. 6 km/sec. For this case, the value of f, is 1. 07. 

s t 

That is, for a particular meteoroid mass, there will be a 7 percent greater meteoroid 
flux on the orbiting spacecraft as compared with the meteoroid flux detected on a sta- 
tionary spacecraft near the Earth. 


Total Meteoroid Velocity Distribution Relative to an 
Orbiting Spacecraft 

In the section entitled Velocity Distribution From Uniform-Velocity Meteoroids, 
it was shown that a velocity distribution develops relative to an orbiting spacecraft from 
meteoroids having uniform velocity V relative to the Earth. This velocity distribu- 
tion is expressed in equation (46). m 

In the section entitled Increase in Meteoroid Flux on an Earth- Orbiting Space- 
craft, the flux of meteoroids on the orbiting spacecraft was found to increase by the 
factor of f(V m ) (eq. (53)) over the meteoroid flux detected on a stationary spacecraft 

at the same position. Thus, the total velocity distribution relative to the spacecraft 
can be obtained by multiplying the fraction of the meteoroids with velocity V m by 

n g (V) and f for all values of V m and then summing the numbers having equal values 

of V. 


However, for most applications, all that is required is to obtain a weighted- 

average velocity. The weighted-average velocity V n (V ) from the distribution of 

n (V) is given in equation (49). Thus, given a distribution of V , the resulting 
s m 

weighted-average velocity relative to the spacecraft can be found from 


/"v,m( V J £ ( V n 

»J dV m 

/ V 

V m) f 

(V )dV 

V m; m 


(57) 


The integral in the denominator is f^, as defined in equation (55) and evaluated in 
equation (56). Thus, if equations (55), (53), and (49) are combined with equation (57) 


V = 


2(n+3)V s f 


;/ 


"V, m( V m)V 2 L( V m + V s f ' ( V m ' dV 


m+3 


\n+3 


m 


(58) 
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Of course, equation (58) can be numerically integrated for any value of n, but it is 

-!/n 

simpler to assume integer values and produce a plot of V as a function of n for 

various values of V . Thus, if n = -2, -1, 1, 2, 3, and 4, the following equations are 
obtained. 




(v 1 ) =f t ( 


m 


-1 


V + v 2 v _1 

fr _ m s m 
» — £ 


(p) 1 ' 2 = ( V m 2 1 2V s 2 +S VV WZ 

f t 


( 7) 1/3 = ( 

f J. 


a/3 


r- | ;V/4 

kV7 


V m 4+5V S 2v m 2 * 3V s 4+ 7 V S V 2 
T 


a /4 


(59) 

(60) 

(61) 

(62) 

(63) 

(64) 


where f^ is given by equation (56). 

If the velocity distribution in reference 3 is used, values for V m n are obtained 

from figure 4. If values for V of 0, 5, 8, and 10 km/sec are assumed, the family 

s 

of curves shown in figure 6 results. Note that the curve for V = 0 in figure 6 is the 
same as the curve in figure 4. The value of n is determined from equation (30). 
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10 12 3 

Power to which velocities are averaged, n 

Figure 6. - Weighted-average meteoroid 
velocity relative to the spacecraft as a 
function of weighting for various space - 
craft velocities (velocity distribution 
from ref. 3). 


METEOROID FLUX AND VELOCITY 
AS FUNCTIONS OF 
DISTANCE FROM THE EARTH 

As a result of Earth gravity and 
shielding, both meteoroid flux and mete- 
oroid velocity distribution vary as func- 
tions of distance from the Earth. In this 
section, these variations are discussed 
relative to the Earth. Velocities relative 
to a moving spacecraft are found in the 
section entitled. Meteoroid Flux and Impact 
Velocity Distribution for an Orbiting 
Spacecraft, after the velocity distribution 
relative to the Earth is known. 

Flux Variation 

Earth shielding factor . - The Earth 
shielding factor is usually calculated by 
taking the fraction of solid angle that the 
Earth subtends and subtracting the frac- 
tion from 1, or 


V = 


1 + cos ip 


(65) 


where 



Figure 7. - Earth shielding factor. 


r 

sin ip = -y (66) 

and r and r are the Earth radius and 
e 

the distance from the center of the Earth, 
respectively (fig. 7). However, for equa- 
tion (66), the assumption is made that 
meteoroids are traveling in a straight line 
between the Earth and r; this, of course, 
is the limiting case as a meteoroid veloc- 
ity increases. 

The effective shielding angle ip e is 

found (through conservation of angular mo- 
mentum (ref. 6)) to be 


sin ip 

e 


V_ r 

oo e 

V r 
m 


(67) 


20 



where V is the velocity of the meteoroid at r , and V is the velocity of the 
00 J e m 

meteoroid at r. Conservation of energy gives V m as a function of V M and r, 


or 



( 68 ) 


Thus, the effective shielding factor becomes 


1 + cos 

\ 2— s (69) 


Since 77 is a function of meteoroid 
6 

velocity, 17 must be integrated over the 

meteoroid velocity distribution in order to 
obtain the total shielding factor. However 
by examining the ratio V e /v for various 

meteoroid velocities (as is done in fig. 8), 
it can be seen that the maximum error in- 
troduced by using rj is only 10 percent, 
at most, and that for most velocity distri- 
butions and distances from the Earth, the 
error is less than 5 percent. Thus, de- 
pending on the accuracy required, equa- 
tions (65) and (66) will usually describe 
the Earth shielding factor accurately. 



Figure 8. - Comparison of Earth 
shielding factors. 


Gravitation decrease factor. - Mete- 
oroid flux will decrease as distance from 
the Earth increases, as a result of the inverse variation of the Earth gravitational field. 
In reference 1, this factor is determined from conservation of angular momentum, and 
the factor can be expressed as 



If equation (70) is combined with equation (68) 


G = 



(71) 


To determine the total decrease in flux resulting from all values of V^, equa- 
tion (71) must be integrated over a distribution of V w , or 


G t = / Gn »( v «) dv « (72 > 


which gives 


G t - 1 - Ve 2 ^ - < 73 > 


!/2 

As an example, when J - 16. 8 km/sec (from fig. 4), 


G, = 1 - 0.44 1 - — 


(74) 


Meteoroid Variation Relative to the Earth 

Because of Earth gravity, the velocity distribution relative to the velocity distribu- 
tion detected on the Earth changes with distance from the Earth. With increasing dis- 
tance from the Earth, meteoroid velocities are slower, but the relative number of 
slower meteoroids also decreases. The net effect is to change some weighted-average 
velocities only slightly, to increase the weighted averages that weight the higher veloc- 
ities, and to decrease the weighted averages that weight the lower velocities. 

If the velocity distribution at the Earth surface, relative to the Earth, is given 
by n oo (V oo ), then the velocity distribution at a distance r from the Earth is given by 



(75) 
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where G is given in equation (71) and is related to V m inequation (68). By 
substituting equations (68) and (71) into equation (75), (V m ) can be determined. 

( n\ l/n 

However, since all that is usually required is determination of (V ) , a determi- 
nation of the new distribution function is not necessary. ' m 

Since 


n f n v m (V m )V m n dV m 

y n _J v, mV m/ m m 

m TnZ (V \dV 

J v , m{ m ) m 

equation (75) is combined with equation (76) to give 

rr~n / G v n ".( v «) dV » 

“ ' ~p ^(V^dV^ 


(76) 


(77) 


Substitution of equations (68), (71), and (72) into equation (77) gives 


V 


n 

m 




-i(n+2)/2 


(78) 


where is given by equation (73). 

Thus, it is necessary to integrate only over the meteoroid velocity distribution 
at the Earth surface. Equation (78) can be integrated directly for values of n = -2, 2, 
and 4. Thus 
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As an example of the application of 
equations (79), (80), and (81), by using 

( iiW n 

values for given in figure 4, the 


three average velocities in equations (79), 
(80), and (81) are found for various dis- 
tances from the Earth. These velocities 
are then plotted as a function of n, and a 
smooth curve is fitted to each set of three 
points. The results are shown in figure 9. 
Three points may seem insufficient to de- 


termine a curve accurately, but when 


intermediate values of 



are 


determined by numerically integrating 
equation (78), the resulting curve proves 
sufficiently well behaved so that three 
points may be considered to be adequate. 


To determine the weighted-average 
velocity of a meteoroid relative to an or- 
biting spacecraft at some distance from 
the Earth, equations (59), (62), and (64) 

can be used, where the values of V n 
’ m 


are found at a distance r by using equa- 
tions (79), (80), and (81). The value of 



for other values of n can be 


— \l/n 

found by plotting \ V/ as a function 


of n. 



Power to which velocities are averaged, n 


Figure 9. - Weighted-average meteoroid 
velocity relative to the Earth as a 
function of weighting for various dis- 
tances from the Earth (velocity dis- 
tribution from ref. 3). 
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CRITERIA 


For the calculations made in this report, the assumption was made that the mete- 
oroid environment at the edge of the Earth atmosphere and relative to the Earth is 
known. The meteoroid flux at the edge of the Earth atmosphere can be expressed as 


F = Cm 


-a 


(17) 


where F is the number of meteoroids of mass m and larger impacting the atmosphere 
per unit area per unit time; C and a are constants. The velocity distribution rela- 
tive to the Earth at the edge of the Earth atmosphere is expressed as n 0O (V 0O ), and the 


meteoroid mass-density distribution is given by n (p). The distribution of impact 

P 

angles into the Earth atmosphere, averaged over the entire surface of the Earth, is 
random, or 


n ( 0 ) = 2 sin 0 cos 0 

0 


( 8 ) 


where 0 is the angle between the meteoroid negative velocity vector and the normal 
vector to the point of impact on the Earth atmosphere. Each of these three distribu- 
tions is obtained by comparing meteoroids of the same mass. 

By using the following method, these distributions may be transformed to a mete- 
oroid flux as a function of some other parameter (e. g. , penetration thickness, energy, 
or momentum) on an orbiting spacecraft at some distance r from the Earth. The 
assumption is made that this other parameter t can be expressed as 


t = Bm^V^p^cos € 0 


(15) 


where m is the impacting mass, V is the relative velocity between the meteoroid 
and impacted spacecraft, p is the mass density of the impacting mass, 0 is the 
angle between the negative velocity vector and the normal to the spacecraft surface at 
the point of impact, and B, 0, y, 5, and e are constants that are characteristic 
of the particular parameter. 

The flux of meteoroids having parameter t and greater is given by 


F = 


G Tjt CB a ^t~ a ^V ya ^ 
t t a 


a a 


(82) 
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where G t is the gravitational decrease factor, 17 is the Earth shielding factor, f^ 
is the flux increase factor on Earth- orbiting spacecraft, V is the weighted-average 

3 

velocity of meteoroids relative to the spacecraft, p is the weighted-average meteor- 
oid mass density, and 6 is the weighted-average meteoroid impact angle. 

Weighted-Average Impact Angle 

The weighted-average impact angle on a randomly tumbling surface is given by 


,/3/ea 


cosV l'^l) 


( 34 ) 


Weighted-Average Meteoroid Mass Density 

The weighted-average mass density is given by 




(3 /5 a 


( 31 ) 


Weighted-Average Impact Velocity 

The weighted-average velocity of meteoroids relative to an orbiting spacecraft is 


=(v n ) 




( 83 ) 


where 


n - 


ya 

T 


( 35 ) 


/— \l/n 

The process of determining \V / from a distribution of velocity relative to the 
Earth is described in the criteria section entitled Determination of the Weighted- 
Average Meteoroid Velocity Relative to an Orbiting Spacecraft at Some Distance From 
the Earth. 
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Gravitation and Earth Shielding Decrease Factors 

The meteoroid mass flux at a distance r from the Earth is decreased by the 
factors 


I 


and 

i 


1 - V V - 2 

0 O0 



(73) 


1 + cos rjy 
2 


(65) 


where 


sin i/' = — 
r 


( 66 ) 


and where V g is the escape velocity from the surface of the Earth, r g is the radius 

of the Earth, r is the distance from the center of the Earth, and ^ is the half angle 
subtended by the Earth from a distance r (fig. 7). 


Flux Increase Factor for an Orbiting Spacecraft 

The meteoroid mass flux on an orbiting spacecraft is increased over the flux on 
a hypothetical stationary spacecraft at the same position by the factor 


f t =1 + 


v 2 v ‘ 

s m 

3 — 


(56) 


where V g is the velocity of the spacecraft relative to the Earth, and V m is the 
velocity relative to the Earth of a meteoroid at a distance r. 
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The weighted average V m 


is related to V 

OO 


by 


V 


-2 

m 



( 84 ) 


Determination of Weighted-Average Meteoroid 
Velocity Relative to an Orbiting Spacecraft 
at Some Distance From the Earth 

Weighted-average meteoroid velocity relative to an orbiting spacecraft . - If the 
meteoroid velocity distribution at a distance r from the Earth and relative to the 
Earth is known to be n^. m (V m ), then the representative velocity of the meteoroids 

relative to the spacecraft is 


V n = 


2(n 






n V, m( V m, 


-2 


)V ( 


m 




( 58 ) 


or, for particular values of n 



( 59 ) 

( 60 ) 

( 61 ) 

( 62 ) 
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( 63 ) 




r 


a/3 



4 9 

V * + 5V Q *V m 
m s m 


W^VsV 2 


a/4 


a 


(64) 


where V is the meteoroid velocity at distance r from the Earth and relative to the 
m 

Earth, V is the velocity of the spacecraft relative to the Earth, and f, is the flux 

S l 

increase factor on an orbiting spacecraft (eq. (56)). 


Usually, it is not necessary to use equation (58) to obtain \V n / for interme- 

diate values of n. The plot of \V / as a function of n is sufficiently well behaved 
that the curve can be accurately interpolated by using equations (59) to (64). An ex- 
ample is shown in figure 6, where it is assumed that the spacecraft is just above the 
Earth surface so that n^ m (V m ) = n oo (V oo ). For this example, the velocity distribu- 
tion at the Earth surface n (V„J is taken from reference 3. If n T7 . (V ) is not 
known, but n 00 (V 00 ) is known, the following section will describe the process to deter- 
mine V n from n (V ). 

oo' oo ' 


Weighted-average meteoroid velocity at some distance from the Earth . - The 
weighted-average meteoroid velocity at a distance r from the Earth and relative to 
the Earth is given by 



•Al' 


f -2 

_ _ 2 

TT 2 

L 


r« 

V 

OO 

V e 

l 1 

-t)\ 


1 (n+2)/2 


”„(v„ 




(’8) 


or, for particular values of n 



(79) 
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where is the meteoroid velocity relative to the Earth at a distance r = r from 

the Earth. (Values of are found from n 00 (V oo ), V e is the escape velocity from 

the surface of the Earth, r is the distance of the spacecraft from the center of the 
Earth, r g is the radius of the Earth, and G t is the gravitational decrease factor 

(eq. (73).) 


as a function of n. 


Equations (79) to (81) are often sufficient to plot 1 V 


Thus, intermediate values of the weighted-average meteoroid velocity can be obtained 
without using equation (78). As an example, figure 9 is a plot in which the meteoroid 
velocity distribution from reference 3 was used. 


A summary of the procedure to find V V / 


from n (V 1 is as follows. 

on ' on ' 


-2 2 4 

1. Numerically obtain the averages , V x , and from n 0O (V 0O ). 


2. Use equations (79) to (81) (and, if necessary, eq. (78)) to plot f V 1 
a function of n for the particular spacecraft distance from the Earth. ' m 


3. Substitute the values obtained for V , 

plot \ V / as a function of n for the particular spacecraft velocity. 


( V n ) 


into equations (59) to (64), and 


4. Read from this plot the weighted-average velocity that corresponds to the 
appropriate value of n (given in eq. (35)). 
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CONCLUSION 


Errors are generated by using average meteoroid velocities and mass densities 
in transforming a meteoroid flux as a function of mass into a flux as a function of pen- 
etration thickness at some distance from the Earth. While in some cases the error 
may be tolerable, in other cases it would be critical. In any event, the use of a prop- 
erly weighted-average velocity is no more difficult than using the average velocity, once 
the representative or weighted-average velocities have been determined for the par- 
ticular spacecraft distance from the Earth and spacecraft velocity relative to the Earth. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, May 27, 1971 
124-12-10-00-72 
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APPENDIX 

ERROR IN METEOROI D VELOCITY AND FLUX DUE TO A 
DIRECTIONAL ENVIRONMENT 


All meteoroid-velocity calculations in this report assume that the spacecraft is 
moving through a meteoroid flux that is random in direction. It is shown that even a 
directional flux becomes random when averaged over the entire Earth surface. How- 
ever, a spacecraft will rarely spend equal time over all areas of the Earth. Instead, 
a spacecraft will circle the Earth in a plane, and its path over the Earth will be a cir- 
cle around the Earth. Although any directional effects will tend to average out along 
the spacecraft trajectory, it may be useful to know what errors are introduced when a 
random directional f lux is adopted. 

The effects of a directional meteoroid stream during one particular spacecraft 
orbit can be quickly calculated. Assume that all meteoroids impact the Earth from a 
single direction in space and that a spacecraft is orbiting the Earth such that the plane 
of its orbit is perpendicular to the direction of the meteoroid stream. Meteoroids 
would impact the spacecraft with velocity V given by 


V = (v 2 + V 2 ) ^ (A1 ^ 

\ m s / 


where is the meteoroid velocity relative to the Earth, and V g is the spacecraft 
velocity relative to the Earth. 

The ratio of the flux on the moving spacecraft to the flux on a spacecraft station- 
ary relative to the Earth is then given by 


f = 




<A2) 


(See the section entitled Velocity Distribution From Uniform- Velocity Meteoroids. ) If 
equation (A2) is integrated over the distribution of V m (eq. (55) from the text) 


‘f/' n ( V m ) dV ffl 


(A3) 
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and the weighted-average meteoroid impact velocity on the spacecraft is (by using 
eq. (57) from the text) 



V 




(n+l)/2 


dV 


m 


(A4) 


If equation (A3) is integrated for n = -1, 


1, and 3 



(A5) 




(A6) 


— / V 3 + 2V h + V 4 V -1 

m s s m 


, 1/3 


(A7) 


Equations (A5) and (A6) are identical to equations (60) and (61), respectively, from the 
text, while equation (A7) is almost identical to equation (63) from the text. Thus, for 
these particular cases, the value of f t represents the principal differences in these 

weighted-average meteoroid velocities. 

If equation (A3) is numerically integrated over the meteoroid velocity distribu- 
tion given by reference 3, and if a value for V is assumed to be 7. 6 km/sec, then a 

s 

value of 1. 096 for f^ is found (compared with a value of 1. 068 for f^ found in the 

section entitled Increase in Meteoroid Flux Relative to an Earth- Orbiting Spacecraft). 
Thus, the error in meteoroid velocity and flux is only 2. 6 percent, even in the unlikely 
circumstance that such a highly directional environment is encountered. The error 
becomes less as V becomes smaller and as the high degree of directionality is 
reduced. s 
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